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Abstract 

The decay B° — > DK*° and the charge conjugate mode are studied using l.Ofb -1 
of pp collision data collected by the LHCb experiment at ^fs = 7TeV in 2011. The 
CP asymmetry between the B° — > DK*° and the B° — > DK*° decay rates, with the 
neutral D meson in the CP-even final state K + K~ , is found to be 

A$ K = -0.45 ±0.23 ±0.02, 

where the first uncertainty is statistical and the second is systematic. In addition, 
favoured B° — > DK*° decays are reconstructed with the D meson in the non-CP 
eigenstate K + ir~ . The ratio of the £>-flavour averaged decay rates in D decays to 
CP and non-CP eigenstates is measured to be 

1l$ K = 1.36 ±ojz± 0.07, 

where the ratio of the branching fractions of D° — > K~tt + to D° — > K + K~ decays 
is included as multiplicative factor. The CP asymmetries measured with two control 
channels, the favoured B° -)■ DK*° decay with D ->■ K + ir~ and the B® ->■ DK*° 
decay with D — > K + K~, are also reported. 
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1 Introduction 



Direct CP violation can arise in B° — > DK*° decays^ from the interference between 
two colour-suppressed transitions: b — > c (Cabibbo favoured) and b — >■ u (Cabibbo 
suppressed). The corresponding Feynman diagrams are shown in Fig. [T| interference 
occurs if the D° and D° mesons decay to a common final state. The magnitude of the CP- 
violating asymmetry that arises from this interference is related to the value of the phase 
7 = arg [— (V u dV* b ) / (y c dV* b )}, the least- well determined angle of the Unitarity Triangle. A 
method to determine 7 from hadronic -B-decay rates was originally proposed by Gronau, 
London and Wyler (GLW) in Ref. |1| for various charged and neutral B — > DK decay 
modes and can be applied to the decay mode B° — > DK*°. In this mode, the charge of 
the kaon from the K*° — > K + ir~ decay unambiguously tags the flavour of the decaying B 
meson j2], hence no time-dependent tagged analysis is required. 

The use of these neutral B decays is particularly interesting because the magnitude 
of the ratio of the suppressed over the favoured amplitude, which controls the size of 
the interference, is expected to be relatively large (naively a factor three larger than the 
analogous ratio for B + — > DK + decays), hence the system can exhibit large CP- violating 
effects, depending on the D decay. Among the modes used in the GLW method, which 
are studied in this paper, large CP asymmetries can be expected when the D meson 
is reconstructed in a CP eigenstate. Contributions from B° decays to the non-resonant 
DK + ix~ final state, which can pollute the DK*° reconstructed signal candidates due to 
the large natural width of the K*°, can be treated in a model- independent way, as shown 
in Ref. |3]. Studies with simulated events have shown that the B° — > DK*° mode is one of 
the most promising channels to provide a precise measurement of 7 at LHCb [4J. Results 
with this channel will therefore complement those from B + — > DK + ', which have recently 




Figure 1: Feynman diagrams for (a) B° -> D°K*° and (b) B° -> D ^* . 



1 Here and in the following, D represents a neutral meson that is an admixture of D° and D°. Inclusion 
of charge conjugate modes is implied unless specified otherwise. 
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been used by LHCb to constrain 7 [HJ. 

This paper presents the measurement of the B° — B° partial width asymmetry using 
D decays into the CP eigenstate K + K~ , 

a kk = r(£° -» D [k +k-]K*) ~ r(B° -» D[K + K-]K*°) 
d " F(B° -> D [K+K - } K* ) + r(£?° D^+jf-]^) ' 

together with the measurement of the ratio of the average of the B° and B° partial widths 
with D — y K + K~ , to the average partial width with D — y K + n~ (where the sign of the 
kaon charge from the D decay is the same as that of the kaon from the K*° decay), 

u kk = r(£° -> D [k+ k-]K*°) + T(B° -> D [k+ k-]K*°) 
d r(i?o *<>) + r(5° -> D [Jf+ff -]A-*o) ' 

These quantities can be used together with other inputs to determine the value of 7. Note 
that the suppressed decay mode B° — ► D^-^+jiT* , where the sign of the kaon charge 
from the D decay is opposite to that of the kaon from the K*° decay, is not included in 
this analysis. This decay mode can exhibit large CP-violating effects and can be studied 
with a larger dataset. The measured asymmetry in the favoured decay B° — y D[K+n-]K*°, 

^fav = r(i? -> Dy-^K*) - T(B° -» Di K + 7T -]K*°) 
d " T(B° -> D [K -^K*°) + r(i?o D[ K + W -]K*°) 

is a useful cross-check since it is expected to be compatible with zero given the size of the 
current dataset. 

In pp collisions, B° s mesons are produced and can decay to the same final state, 
B° s — y DK*° |6j. In these B° s decay modes, the interference between the two contributing 
amplitudes is expected to be small, since the relative magnitude of the suppressed to the 
favoured amplitude is small compared to the B° case. Therefore, these modes are valuable 
control channels, and the asymmetry 

a kk = r(g -» D [k+ k-]K*°) - T(B° S -> D[ K + K -]K*°) 
s ~ r(5° -> D^+k^K* ) + r(B° D[ K + K -]K*°) ' 

similar to that defined in Eq. (lip, is also obtained in this analysis. Since the favoured 
(suppressed) B° s (B°) decay gives kaons with opposite charges from D and K*° decays, 
Ag is not used as a control measurement in the analysis, to avoid biasing a potential 
future measurement of „4!, up . 



2 The LHCb detector, dataset and event selection 

The study reported here is based on a data sample collected at the Large Hadron Collider 
(LHC) with the LHCb detector at a centre-of-mass energy of 7TeV during the year 2011, 
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corresponding to an integrated luminosity of 1.0 fb~ . The LHCb detector [7] is a single- 
arm forward spectrometer covering the pseudorapidity range 2 < rj < 5, designed for the 
study of particles containing b or c quarks. The detector includes a high precision tracking 
system consisting of a silicon-strip vertex detector surrounding the pp interaction region, 
a large-area silicon-strip detector located upstream of a dipole magnet with a bending 
power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes 
placed downstream. The combined tracking system has a momentum resolution Ap/p that 
varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and an impact parameter resolution of 
20 \mi for tracks with high transverse momentum (pt). Charged hadrons are identified 
using two ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are 
identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers. The 
trigger (8) consists of a hardware stage, based on information from the calorimeter and 
muon systems, followed by a software stage which applies a full event reconstruction. 

This analysis uses events selected by the hardware level trigger either when one of the 
charged particles of the signal decay gives a large enough energy deposit in the calorimeter 
system (hadron trigger), or when one of the particles in the event, not coming from the 
signal decay, fulfills the trigger requirements (i.e. mainly events triggered by one particle 
coming from the decay of the other B in the event). The software trigger requires a two-, 
three- or four-track secondary vertex with a high scalar sum of the of the tracks and 
a significant displacement from the primary pp interaction vertices (PVs). At least one 
track should have > 1.7GeV/c and an impact parameter (IP) x 2 with respect to the 
PV greater than 16. The IP \ 2 is defined as the difference between the x 2 of the PV 
reconstructed with and without the considered track. A multivariate algorithm is used for 
the identification of secondary vertices consistent with the decay of a b hadron. 

Candidates are selected from combinations of charged particles. D mesons are recon- 
structed in the decay modes D — > K + tt~ and K + K~ . The pr of the daughters is required 
to be larger than 400MeV/c. Particle identification (PID) is used to distinguish between 
charged pions and kaons. The difference between the log-likelihoods of the kaon and pion 
hypotheses (DLL^) is required to be larger than for kaons and smaller than 4 for pions. 
This aids the reduction of cross-feed between the signal D decay modes to a negligible level. 
A fit is applied to the two-track vertex, requiring that the corresponding x 2 P er degree of 
freedom is less than 5. In order to separate D mesons coming from a B decay from those 
produced at the PV, the D candidates are required to have an IP x 2 greater than 4 with 
respect to any PV. To suppress background from B decays without an intermediate D 
meson (B° — > K*°K + K~ for example), for which all four charged hadrons are produced 
at the .B-decay vertex, a condition on the D flight distance with respect to the B vertex is 
applied, requiring that it is larger than by at least 2.5 standard deviations. Finally, D 
candidates with an invariant mass within ±20 MeV/ c 2 of the nominal D° mass are retained. 

K*° mesons are reconstructed in the mode K*° — > K + tt~. The pt of the K + and 
7T~ mesons must be larger than 300MeV/c. PID is also used, requiring that DLLxtt is 
larger than 3 for the kaon and lower than 3 for the pion, reducing the cross-feed from 
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B° Dp to a manageable level and rejecting non- resonant B° -»■ DK + K~ |9|. Possible 
contamination from protons in the kaon sample, e.g. from A® — )■ Dp7t~ decays, is reduced 
by removing kaon candidates with a difference between the log-likelihoods of the proton 
and kaon hypotheses (DLL p k) of less than 10. The IP \ 2 °f the K*° mesons must be 
larger than 25, to select those coming from a B decay, and their invariant mass within 

±50 MeV/ c 2 of the nominal mass. 

(— ) 

B( s j meson candidates are formed by combining D and K candidates selected with 
the above requirements. A fit to a common vertex is performed, keeping only combinations 
with x 2 P er degree of freedom lower than 4, and a kinematic fit is performed to constrain 



the invariant mass of the reconstructed D to the nominal D° mass 10 . Since B mesons 
are produced at the PV, only candidates with IP \ 2 lower than 9 are retained. In case 
several PVs are reconstructed, the one for which the 5-candidate IP x 2 is the smallest 
is taken as reference. Additionally, the momentum of the reconstructed B candidate is 
required to point back to the PV, by requiring that the angle between the B momentum 
direction and its direction of flight from the PV is smaller than 10 mrad. Furthermore, the 
sum of the square roots of the IP \ 2 °f the four charged particles must be larger than 32. 
The absolute value of the cosine of the K*° helicity angle is required to be larger than 0.4. 
This angle is defined as the angle between the kaon-daughter momentum direction in the 
K*° rest frame, and the K*° direction in the B rest frame. 

Specific peaking backgrounds from BP*. — > D^Ji decays, where h is a tt or a K meson, 
are eliminated by vetoing candidates for which the invariant mass of K + K~ it + (K~ n + n + 
and K + K~7t + ) is within ±15 MeV/ c 2 of the nominal mass of a {D + ) meson. 

Where possible, data-driven methods are used to determine selection efficiencies and 
invariant mass distribution shapes. Otherwise, they are determined from fully simulated 
events. The pp collisions are generated using Pythia 6.4 111] with a specific LHCb 



configuration 12 where, in particular, decays of hadronic particles are described by 



EvtGen 13 . The interaction of the generated particles with the detector and its 



response are implemented using the Geant4 toolkit 14 as described in Ref. 15 



3 Determination of signal yields 

The numbers of reconstructed signal B° and B® candidates are determined from an 
unbinned maximum likelihood fit to their mass distributions. Candidates are split into four 
categories, which are fitted simultaneously: D(K + R-)K* , D(K + R-)K* , D(K + ir-)K*°, 
and D(K~7T + )K*°. The mass distribution of each category is fitted with a sum of 
probability density functions (PDF) modelling the different contributing components: 

1. the B° and B° s signals are described by double Gaussian functions; 

2. the combinatorial background is described by an exponential function; 

3. the cross-feed from B° — > Dp decays, where one pion from the p° — > tt + 7t~ decay is 
misidentified as a kaon, is described by a non-parametric PDF (16) determined from 
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fully simulated and selected events; 

4. the partially reconstructed B° — > D*K*° and B° s — > D*K*° decays, where the D* is 
a D*° or a D*° and the ir° or photon from the D* decay is not reconstructed, are 
modelled by a non-parametric PDF determined from fully simulated and selected 
events. 

There are 23 free parameters in the fit. These include the B° PDF peak position, 
the core Gaussian resolution for the B° and the B® and the slope of the combinatorial 
background, all of which are common to the four fit categories. The remaining free 
parameters are yields for each fit component within each category. Yields for B^ and B9\ 
are constrained to be identical for the background components where CP violation effects 
can be excluded or are expected to be compatible with zero with the current data sample 
size. 

A separate fit to B° — > D(K + ir~)p candidates in the same data sample is performed. 
The yield of such candidates and the probability to reconstruct them as B° — > D(K + tt~)K*° 
is used to constrain the number of cross-feed events in the D(K + tc~)K*° category. The 
number of cross-feed candidates from B° — > D(K + K~)p° in the D(K + K~)K*° category 
is derived from the D(K + tt~)K*° category using the relative D branching fractions and 
B selection efficiencies. As no flavour asymmetry is expected for this background, the 
numbers of cross-feed events in the DK*° categories are constrained to be identical to 
those of the corresponding DK*° categories. 

The partially reconstructed background component accumulates at masses lower than 
the nominal B° mass. Its shape depends on the unknown fraction of transverse polarisation 

in the B?^ — > D*K*° decays. In order to model the b\\ — > D*K*° contribution, a PDF is 
built from a linear combination of three non-parametric functions corresponding to the three 

orthogonal helicity eigenstates. The functions are derived from simulated B? s \ — > D*K*° 
events reconstructed as B° — > DK*°. Each function corresponds to the weighted sum of 
the D* — > D'y and D* — > Dtt° contributions for a defined helicity eigenstate, where the 
weights take into account the relative D* decay branching fractions and the corresponding 
reconstruction efficiencies. 

The invariant mass distributions together with the function resulting from the fit are 
shown in Fig. [2] Note that the decay B° s — > D(K + 7r~)K*° is not observed since the charge 
combination of the kaons in the final state corresponds to the suppressed decay. The signal 
yield in each category is summarized in Table [!} The significance of the B° — > DK*° 
signal for D — > K + K~ decays, summing B° and B° and including both statistical and 
systematic uncertainties, is found to be equal to 5.1 cr, by comparing the maximum of the 
likelihood of the nominal fit and the maximum with the yield of the B° — > D(K + K~)K*° 
category set to zero. 

The yields determined from the simultaneous mass fit are corrected for selection 
efficiencies in order to evaluate the asymmetries and ratios described in the introduction. 
The selection efficiencies account for the geometrical acceptance of the detector, the 
reconstruction, the PID, and the trigger efficiencies. All efficiencies are computed from 
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fully simulated events, except for the PID and trigger efficiencies, which are obtained 
directly from data using clean calibration samples of D° — > K~n + from D* + decays. 



Table 1: Signal yields with their statistical uncertainties. 
Category Signal yield Category Signal yield 



B° -> D 
B° D 



[K + TT-] 



]K *o 


21 




23 


K*° 


108 



- 5 

- 12 
11 



B° -> D^+k^K* 
B° s -> D [K+K - } K*° 
B° -> D [K - n+] K*° 



8±4 
24 i« 
94 ±11 
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Figure 2: Invariant mass distributions of (a) Di K + K -iK*°, (b) Di K + k --iK*°, (c) D[ K - n +^K*° and 
(d) Drj(+ w -iK* candidates. The DK*° distributions correspond to B° and B® decays whereas 
the DK*° distributions correspond to B° and B® decays. The fit functions are superimposed; 
the different B decays and combinatorial background components are detailed in the legends. 
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4 Systematic uncertainties 



Several sources of systematic uncertainty are considered, affecting either the determination 
of the signal yields or the computation of the efficiencies. They are summarized in Table [2j 
In order to take into account the measured difference in the production rate between B° 
and B°, the B° yields are multiplied by a correction factor, 

d _ 1 - KA prod 
"prod - 1 a , l°J 



where v4 pro( j = 0.010 ± 0.013 17 is the asymmetry between B° and B° at production in 



n . . , . , , , , c \ ;+°° e~ Tt cos(Ami) e(B°-s>_D_ft:* ,t) dt n . , 

collisions, and k is a decay-dependent factor, k = — — +oc g _ rt ^o^^^o ^ df > which 

takes into account dilution effects due to the B°—B° oscillation frequency, Am, and includes 
the acceptance as a function of the decay time for the reconstructed decay, e(B° — y DK*°, t). 
The value of k is found to be 0.46±0.01 using fully simulated events and PID efficiencies from 
calibration samples. The uncertainty on ap rod is propagated to the measured observables 
to estimate the systematic uncertainty from the production asymmetry. Owing to the 
large B° s oscillation frequency, the potential production asymmetry of B® mesons does not 
significantly affect the measurement presented here and is neglected. 

The PID calibration introduces a systematic uncertainty on the calculated PID effi- 
ciencies, which propagates to the final results. All PID correction factors are compatible 
with unity within their uncertainties which are of the order of 1%. 

The systematic uncertainty associated to the trigger is estimated by varying in the 
simulation the fraction of events triggered by the hadron trigger with respect to the fraction 
of events triggered by the other &-hadron in the event. Other selection efficiencies cancel in 
the ratio of yields, except for the efficiencies of the pr cuts on the D daughters, which are 
different between different D decay modes. TZ^ K has to be corrected by a multiplicative 
factor 0.94 ± 0.04, where the statistical uncertainty on the correction, which arises from 
finite simulated sample size, is assigned as systematic uncertainty due to the relative 
selection efficiencies. 

The fit procedure is validated with simulated experiments. A bias of statistical nature, 
owing to the small number of events in the B° — > D(K + K~)K*° channel, is found to be 
5% for B° and 8% for B°. The signal yields are corrected for this bias before computing the 
asymmetries and ratios. A systematic uncertainty equal to half the size of the correction 
has been assigned. 

Simulated experiments are also used to determine the systematic uncertainties due to 
the low-mass background, the B° — y Dp cross-feed, and the signal shape. Samples are 
generated with different values of the polarisation parameters, the cross- feed fraction and 
the fixed signal parameters. The corresponding systematic uncertainty is estimated from 
the bias in the results obtained by performing the fit described in the previous section to 
these samples. 
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Table 2: Summary of the absolute systematic uncertainties on the measured observables. 



Source 


AKK 


^d 




IK -d 


Production asymmetry 


0.005 


0.006 




0.003 


PID efficiency 


0.004 


0.008 


0.005 


0.014 


Trigger efficiency 


0.004 


0.001 


0.005 


0.022 


Selection efficiency 








0.040 


Bias correction 


0.004 




0.001 


0.013 


Low-mass background 


0.017 


0.001 


0.004 


0.042 


B° — > Dp cross- feed 


0.001 




0.002 


0.008 


Signal description 


0.001 


0.001 


0.001 


0.005 


D branching fractions 








0.022 


Total 


0.019 


0.010 


0.008 


0.069 



5 Results and summary 

This paper reports the analysis of B° — > DK*° decays using 1.0 fb _1 of pp collision data. 
Potential contributions to the decay amplitudes from the non-resonant B° — > DK + tt~ 
mode are reduced by requiring that the K*° reconstructed mass is within ±50MeV/c 2 of 
the nominal mass and the absolute value of the cosine of the K*° helicity angle is greater 
than 0.4. The results for the CP-violating observables are 

A% K = -0.45 ± 0.23 (stat) ± 0.02 (syst), 
Af v =-0.08 ±0.08 (stat) ±0.01 (syst), 
A* K = 0.04 ± 0.16 (stat) ± 0.01 (syst), 
K% K = 1.36 t oil (stat) ± 0.07 (syst). 

The value of TZ^ K takes into account the ratio of the branching fractions of D° — > K + K~ 
to D° — > K~n + decays 10 . The correlation between A^ K and TZ^ K is equal to 0.16 and 



the correlations between the other observables are negligible. 

These are the first measurements of CP asymmetries in B° and B® to DK*° decays 
with the neutral D meson decaying into a CP-even final state. With more data, improved 
measurements of these and other quantities in B° — > DK*° decays will result in important 
constraints on the angle 7 of the Unitarity Triangle. 
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